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1. Introduction 

The question of coexistence of magnetism and superconductivity appears very often 
in correlated electron systems. In this context, both the spin-singlet and the spin- 
triplet paired states should be considered. A general motivation for considering here 
the spin-triplet pairing is provided by the discoveries of superconductivity in Sr2Ru04 
[H [2], UGe2 [a m, URhCe [5], Ulr |6], and UCoGe [3 El E]. In the last four 
compounds, superconductivity indeed coexists with ferromagnetism. Moreover, for 
both, the spin-singlet high-temperature superconductors and the heavy-fermion systems, 
the antiferromagnetism and the superconductivity can have the same origin. Hence, 
it is natural to ask whether ferromagnetism and spin-triplet superconductivity also 
have the same origin in the itinerant uranium ferromagnets. A related and a very 
nontrivial question is concerned with the coexistence of antiferromagnetism with triplet 
superconducting state as in UNigAlg [lOl [IH [12] and UPtg ^ [H]. 

It has been argued earlier [T5l \T7\ [T8| [T9] that for the case of indistinguishable 
fermions, the intra-atomic Hund's rule exchange can lead in a natural manner to the 
coexistence of spin-triplet superconductivity with magnetic ordering - ferromagnetism 
or antiferromagnetism in the simplest situations. This idea has been elaborated 
subsequently by us [20l [211 [22] by means of the combined Hartree-Fock(HF)-Bardeen- 
Cooper-Shrieffer(BCS) approach. In particular, the phase diagrams have been 
determined which contain regions of stability of the pure superconducting phase of 
type A (i.e., the equal-spin-paired phase), as well as superconductivity coexisting with 
either ferromagnetism or antiferromagnetism. 

The HF approximation, as a rule, overestimates the stability of phases with a 
broken symmetry. Therefore, in this work, we apply the Gutzwiller approximation for 
the same selection of phases in order to examine explicitly the effects of interelectronic 
correlations. The extension of the Gutzwiller method to the multi-band case |23l \2M [25] 
provides us with the so-called renormalization factors for our degenerate two-band 
models. With these factors we construct an effective Hamiltonian by means of the 
statistically consistent Gutzwiller approximation, SGA, in which additional constraints 
are added to the standard Gutzwiller approximation (GA) and with the incorporation 
of which the single-particle state has been determined (see [261 [271 [28l [29] for exemplary 
applications of the SGA method). The detailed phase diagram and the corresponding 
order parameters are determined as functions of the microscopic parameters such as 
the band filling, n, the Hund's rule exchange integral, J, and the Hubbard interaction 
parameters, U and U'. The obtained results are compared with those coming out from 
the Hartree-Fock approximation. In this manner, the paper extends the discussion of 
itinerant magnetism within the canonical (extended Hubbard) model, appropriate for 
this purpose, to the analysis of pure and coexisting superconducting-magnetic states 
within a single unified approach. Additionally, at the end, we dwell briefly on the 
applicability of our original concepts to more realistic systems. It should be noted that 
theoretical investigations regarding the spin-triplet pairing have been performed recently 
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also for other systems [201 ED El 1331 [Ml E5] • 

The paper is composed as follows. In Sec. II we provide the principal aspects of 
real-space spin-triplet pairing induced by the Hund's rule coupling, and introduce the 
band-renormalization factors for our two-band model. Furthermore, in subsections A 
and B of Sec. II we explain how the effective Hamiltonian is constructed, according to 
the statistically consistent Gutzwiller approximation, for all the phases considered in 
this work. In Sec. Ill we discuss the phase diagram, and the principal order parameters 
in the considered phases, whereas Sec. IV contains the concluding remarks and outlook 
concerning future investigations to make the approach applicable to real systems. 

2. Model and method 

We consider the extended orbitally-degenerate Hubbard Hamiltonian, which has the 
form 

+ U^ fiii^nui - J ^ ( Si/- Sw + -huhii> j (1) 

= ^° + H^^ , 

where / = 1,2 label the orbitals and the first term describes electron hopping between 
atomic sites i and j. For / ^ I' this term represents electron hopping with change 
of the orbital (i.e., hybridization in momentum space). The next two terms describe 
the Coulomb interactions between electrons on the same atomic site. However the 
second term contains also the contribution, originating from the exchange interaction 
(J). The last term expresses the Hund's rule i.e., the ferromagnetic exchange between 
electrons localized on the same site, but on different orbitals. This term contributes 
to magnetic coupling and is responsible for the spin-triplet pairing leading to magnetic 
ordering, superconductivity and coexistent magnetic-superconducting phases. In the 
Hamiltonian ([T]), we have disregarded the pair hopping term (^/2) X]z=^;' ^Izt^Li^j''4^«''t 
because it hardly influences the ordered phases which we analyze in this work. In 
our variational method we assume that the correlated state I^&g) of the system can be 
expressed in the following manner 

|*g) = PgI^o) , (2) 

where |\l/o) is the normalized non-correlated state to be determined later and Pq is the 
Gutzwiller correlator selected in the following form 

i i i,r 
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Here, Xjj, are the variational parameters, which are assumed to be reaL In the two-band 
situation the local basis consists of 16 states (see Table H]), which are defined as follows 






(4) 



where 7 = 1, 2, 3, 4 labels the four spin-orbital states (in the la notation: 1 f, 1 J,, 2 f, 2 J,, 
respectively) and \I\ is the number of electrons in the local state \I). In general, an index 
/ can be interpreted as a set in the usual mathematical sense. The creation operators 
in (jll) are placed in ascending order, i.e., 71 < ... < 7|/|. In an analogous manner, one 
can define the product of annihilation operators 



C^,I = ll 



-«7 



Cj7i---Q7|j| ) 



(5) 



76/ 



which are placed in descending order 71 > ... > ^\j\. 



Table 1. The local basis consisting of 16 configurations containing N^ 
electrons, which are enumerated as shown below. 



= 0, 



|0,0) 
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0,1) 
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|i,i) 
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ti,t) 


13 


lt,o) 
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ti,o) 


6 


lt,i) 


10 


i,ti) 


14 


|o,t) 


3 


o,ti) 


7 


li,t) 


11 


ti,i) 


15 


;,o) 


4 


t,t) 


8 


lt,ti) 


12 


ti,ti) 


16 



The operator |/)jj(/'| can be expressed in terms of Cj and C*/ in the following 
manner 

m,i'\i ^ \I)n{I'\ = ClA,i'nU'\i , (6) 

where 



"'iur\i 



n 

7G/U/' 



:i-n 



I'y ) 



(7) 



In the subsequent discussion, we write expectation values with respect to |\E'o) as 

(6)0 = (^o|0|^o) , (8) 

while the expectation values with respect to \^g) will be denoted by 

(*g|0|^g) 



(O) 



G 



(9) 



The most important step within the Gutzwiller approach is to derive the formula for the 
expectation value of the Hamiltonian K = H — fiN with respect to |\E'g)- This can be 
done in the limit of infinite dimensions by a diagrammatic approach |25j which uses the 
variational analog of Feynmann diagrams. By applying this method to the interaction 
part of the Hamiltonian ([1]) , which is completely of intra-site character, one obtains 



{H''')g 



h,l4 



(10) 
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where 



E 






:ii) 



and L is the number of atomic sites. In f ITOl) we have assumed that our system is 
homogeneous. Note that, with the use of Wick's theorem, the purely local expectation 
values {rhi^j^Q can be expressed in terms of the local single-particle density matrix 
elements {cf cf i)q. Here, cf are either creation or annihilation operators. 

The expectation value of the single-particle part in the Hamiltonian ([T]) can be cast 
to the form 

{H°)g= Y1 ^^l]'{^ii^i'i'-^ii^i'i')'\C^in^Cin')o (12) 

«j(«7^i)77'77' 

where we have assumed that the renormalization factors q and q are real numbers and 
^77 — fii _ Moreover, in the equation above we have neglected the part containing 
the inter-site pairing terms (cj,^c^ ,)o and (cj,^Cj,-y/)o as we are going to concentrate 
on the Hund's rule induced intra-site spin-triplet paired states. The inter-site pairing 
amplitudes are much smaller than the intra-site terms, in the considered model. The 
renormalization factors, introduced in (IT2|) . have the form 



Q'77 



^ ^fsgn(7,/)m;j7?c^U7,/'|7+ 5Z ^^§'^(7. ^)"^/',/U7C/',/ 



|7 



imi) "- I' 



I'mn 



(13) 



0(7) 



,(7) 



where m%, = (m/,/')o and nij], = {m\'\,)Q. Here we have introduced the operator 



,{7) 



^t 



The parameters d 



h,l2\^ 



in (IT3I) are defined as 



^*h,i2h= X] fsgn(7,/)A/,,/u7A/,/2 

7(7^/) 



(14) 



(15) 



where we introduced the fermionic sign function 

fsgn(7,/) = (/U7|c;|/). 



(16) 



The renormalization factors q^^ have to be included in flT2|) , when there are nonzero gap 
parameters ((c"c")o 7^ 0) in |\l/o)) which is the case considered here. The form of q.yj is 
as follows 



977 



Y^ ^ fsgn(7, J)m;f J c;^,^jU7l7 + Yl ^^&^(^^^)^W'^li' 



|7 



imn "- I' 



I'mn 



(17) 



The remaining part of {K)g that has to be derived is the expectation value {N)g- Also 
in this case, the diagrammatic evaluation in infinite dimensions gives the proper formula, 

{N)g = Y.^n,,)G , (18) 

27 
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where 







(^n)G 








(19) 


and 




N7 r = 


/ , '^/i,/U7'^/U7,/4 ■ 

7(7^/) 






(20) 


The pairing 


densities in the correlated state that are going 


to be 


useful 


in the subsequent 


discussion can be expressed 


in the following way 









where 



(ci7Ciy)G = ^ Sj^^,ijri\,u ' (21) 

SVm= E A,,,A7u(77'),/4fegn(7,/)fsgn(y,/)fsgn(7',7). (22) 

7(77'^/) 



Using ( TTOj) . ( IT2|) . and ( 1T8l) one can express {K)g in terms of the variational 
parameters A/.//, local and non-local single particle density matrix elements ,{c2yC" i)o, 
{c] Cj^y)o, and the matrix elements of the atomic part of the atomic Hamiltonian 
represented in the local basis (/|i/"*|/'). 

The formula for {K)g, obtained in the way described above, can be written as an 
expectation value of an effective Hamiltonian Kqa, evaluated with respect to |\l'o) 

Kg A = Z^ Z^ ^V (?77'?7'7' ~ ^77^7'7')Cj,7Ci,7' 

ij{i¥=j) 77'77' /„o\ 

-f^Yl ^7^*^ + ^ X] ^h,h{^ 11,14)0, 

where q^ = {hi^)G/ {f^i^)o- There is no guarantee that the condition 

{hi^)G = (^n)o 5 (24) 

is fulfilled. It turns out that it is fulfilled for the paramagnetic and the magnetically 
ordered phases of our two-band system, however it is not for the superconducting phases. 
Physically it is most sensible to fix {n)G instead of (n)o, during the minimization. This 
is the reason why we include the term — /xiV already at the beginning of our derivation in 
{K)g- In this manner the chemical potential /i refers to the initial correlated system, not 
to the effective non-correlated one (for which the chemical potential can be different). 

Having in mind that there are 16 states in the local basis there could be up to 
16 X 16 = 256 variational parameters Xij'. However, for symmetry reasons many of 
these parameters are zero. The finite parameters can be identified by the following rule 

Xi,r + ^ (m7,7,)o + V (/|i^'^*|/') ^ 0; . (25) 
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It should also be noted that, as shown in PS], the variational parameters are not 
independent since they have to obey the constrains 

(^GI.)O = 1, 

/At p2 At \ _ /At At \ ^ ' 

{^ij-^Gli^iy/O = {(^i-yCi'y' ) , 

which are going to be used to fix some of the parameters Xij/. 

The results presented in this work have been obtained for the case of a square lattice 
with the band dispersions 

eik = e2k = ek = 2t(cos (k^^) + cos (ky)) , (27) 

and also 

ei2k = e2ik = /^hCk , (28) 

where Ph G [0, 1]. The orbital degeneracy and spatial homogeneity allow us to write 

{hii)G = {ni2)G = ^g/2, 

/ QZ \ / QZ \ OZ ^ ' 

Wil/G - Wi2/G = ^G ' 

where 

nn = riii'^ + riii^ . 

Similar expressions as in (I29p can be introduced for the expectation values in the non- 
correlated state |\l/o)- 

Before discussing the principal magnetic and/or spin-triplet superconducting 
phases, we introduce first the exact expression of the full exchange operator (the last 
term of our Hamiltonian) via the local spin-triplet pairing operators {A]^, Aim) namely 

y^ ( Si/- Sw + -hiifiii, j = ^ A\^Aim , (31) 

where 

{«llt«l2t "^ = 1 

AxAn "^ = -1 (32) 

We see that those two representations are mathematically equivalent, so the phase 
with Sq = {Sii)G 7^ and that with the corresponding off-diagonal order parameter 
{Aim)G 7^ (or {AI^)g 7^ 0) should be treated on equal footing. 
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2.1. Statistically- consistent Gutzwiller method for superconducting and coexistent 
superconducting-ferromagnetic phases 

In this subsection we will describe the SGA approach as applied to the selected phases 
characterized by the following order parameters 

• Superconducting phase of type Al coexisting with ferromagnetisni (Al+FM): 
5^l„ ^ 0, Af ^ 0, A^, = A^ = 0, 

• Pure type A superconducting phase (A): 



^G\u 



0, Af = A^, ^ 0, A« = 0, 
Pure ferromagnetic phase (FM): 
Shiu 7^ 0, Af = A^, 



A^ 



Paramagnetic phase (NS): 
0, Af = A% 



qz 



AG' 



0, 



0, 



where Sq,^ refers to the uniform magnetic moment and 



A 



G 



(A,^)g, (A^)* = (iL) 



(33) 



are the spin-triplet local gap parameters which are assumed as real here. 

The (correlated) order parameters which have been used above to define the relevant 
phases can also be defined for the non-correlated state I^E'o)- With these, we can 
determine which of the matrix elements (m///)o are equal to zero for the considered 
phases. The assumption (125|) then allows us to choose the non-diagonal variational 
parameters, A/,//, that have to be taken into account during the calculations. We list 
their indexes (J, I') in Table [2l 



Table 2. Nonzero, off-diagonal local variational parameters (A/,// = A/'_/) that are 
used in the calculations for the considered phases. 
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10 
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/' 


16 


15 


14 


13 


12 


16 


16 
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11 


8 


9 



As one can see from Table [21 the off-diagonal variational parameters correspond to 
the creation or annihilation of the Cooper pair in the proper spin-triplet states 1 1 t) 2 t) 
and |1 1,2 l) (phase A). Because in the Al phase only electrons with spin- up are 
paired one can assume that Ai^ie, A2,i5, A3,i4, A8,i6, A8,9 are zero (and their transoposed 
corespondants - Xjj' = Xpj). For the FM and NS unpaired states only Aio,ii and An^io 
are nonzero. They correspond to the two non-diagonal matrix elements of the atomic 
Hamiltonian, {I\H°-^\I'). With the information contained in Table [21 one obtains the 
following relations regarding the band-narrowing renormalization factors 



qia,i'a' ^0^1 = 1' Aa = a', 



(34) 
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where we have again used the '-/ = la notation. Due to the degeneracy of our bands we 
find 

QlaM = C[2a,2cr = Qa , 

(l2cr,lcr = (llcr,2a = Qa , (35) 

qL = qL = € ■ 



Using the equations above we can rewrite the Hamiltonian (1231) in the more explicit 
form, in reciprocal space 

KgA = J^lQaCk - qllA^^Ma + ^ Qa^\,12c\,laCM'a + ^ Xl ^h,lA'>TT'h,u) , (36) 

where the renormalization factors Q^ are defined as 

Q.^ql- € ■ (37) 

Having the formula for KqAi given by (l36l) . one can introduce next the so-called 
statistically- consistent Gutzwiller approximation (SGA). In this method, the mean fields 
(such as the expectation values for magnetization or superconducting gaps) are treated 
as variational mean-field order parameters with respect to which the energy of the system 
is minimized. However, in order to make sure that they coincide with the corresponding 
values calculated self-consistently, additional constraints have to be introduced with the 
help of the Lagrange-multiplier method [26l [27l [28l |29] . This leads to supplementary 
terms in the effective Hamiltonian of the following form 



K^ = Kga- Yl ^- ( Zl ^km - LaA + H.C. 



m=itl 






(3J 



\d ^ ^ V.I 



where the Lagrange multipliers A^, A^, and A„ are introduced to assure that the 
averages (Akm), (•S'kz) and (jivia) calculated either from the corresponding self-consistent 
equations or variationally, coincide with each other |29| . 

Introducing the four-component representation of single-particle operators 

fL = (4ia' c^a' c-kia, c_k2.) , (39) 

we can write down the effective Hamiltonian in the following form 



kcr kcr 7Ti=dil 



+ LYEh,iA^iui4)o 



(40) 
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where Mko- is a 4x4 orthogonal matrix 



10 



M 



kcr 



f ~e.. 


Qa^kU 





A. ^ 


Qa^kU 


eka 


-A. 








-A. 


-CkCT 


— Qa^kU 


iy K 





— Qa^kl2 


-ek<7 ) 



(41) 



Here we introduced A^- and X^ which correspond to the Lagrange parameters Xm=i and 
Xm=-i, respectively. The bare quasiparticle energies e^i^ are defined as 



eka = Qa^k -qa{l^ + A„) - -crXs . 



(42) 



The diagonalization of the matrix ( 14T|) yields the quasiparticle eigen-energies in the 
paired states of the following form 



-E'kio- 


- Y 4,7 + A^ - QcT^kl2 , 


Ek2a 


= \/eL + K + Qaeki2 , 


EkSa 


= -\/eL + A^-g.6ki2 


-E'k4o- 


= -\A~L + A^ + g.eki2 



(43) 



The first two energies correspond to the doubly degenerate spin-split quasiparticle 
excitations in the A phase, whereas the remaining two are their quasihole 
correspondents. 

Even though the Gutzwiller approach was derived for zero temperature, we may 
still construct the grand-potential function Fx (per atomic site) that corresponds to the 
effective Hamiltonian fl^Ol). i.e., 



no 



^^ = -7^ E 1^ [1 + ^"''''"1 + 7 E ^~k. + 2 5^ X^Al + 2XsS^ + (A„ + /i) 



' kla kcr 

^EhM(frihM)o ■ 



m,=itl 



h,lA 



(44) 



The values of the mean fields, the variational parameters, and the Lagrange multipliers 

are found by minimizing the Fx functional, i.e., the necessary conditions for minimum 

are 

riF. dFx dFx , , 

0, ^^ = 0, ^ = 0, (45) 



dA dkv dAi 

where A, Ay, A^ denote collectively the mean fields in the non-correlated state, the 
variational parameters and the Lagrange multipliers respectively. Additionally, the 
chemical potential, /i enters through the relation 

dFx 



dnr 



/i. 



(46) 
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After solving the complete set of equations, one still has to calculate the mean fields in 
the correlated state with the use of their analogs in the non-correlated state and the 
variational parameters using ( IT9|) and ( 12T]) . 

With the SGA method one minimises the variational ground state energy {K)g with 
respect to the variational parameters Xij> and the single-particle states I^E^o)- Note that 
an alternative way for this minimization has been introduced, e.g., in [37|. Beyond the 
ground-state properties of K one is often also interested in the (effective) single-particle 
Hamiltonian (140|) because its eigenvalues are interpreted as quasi-particle excitation 
energies |38] . 

2.2. Statistically- consistent Gutzwiller method for the coexistent 
antiferromagnetic- spin-triplet superconducting phase 

To consider antiferromagnetism in the simplest case, we divide our system into two 
interpenetrating sublattices A and B. In accordance with this division, we define the 
annihilation operators on the sublattices 

J Cii^A ioT ie A , 
1^ Cii^B for 1 G 5 . 

The same holds for the creation operators. Next, the Gutzwiller correlator can be 
expressed in the form 

i{A) i(B) 

where 

^S''^ = E^5TV).(/1. (49) 

1,1' 

If we assume that charge ordering is absent, we have 

WilA/G = Sq^^, {S-iBJG = —Sg\s ) (50) 

{niiA)G = {niiB)G = nG/2 . (51) 

Similar expressions can be obtained for the case of expectation values taken in the 
state I^E'o)- As one can see from ( 149|) . we have introduced separate sets of variational 
parameters {Xfj, and \f ji) for the two sublattices. Fortunately, it does not mean that 
we have twice as many variational parameters as in the preceding subsection. The 
parameters Xf^ are related to the corresponding \f p through 

where the states Ji and I2 have opposite spins to those in the I^ and /4 states, 
respectively. The same division has to be made for the renormalization factors q, q 



Coexistence of spin-triplet superconductivity with magnetism 



12 



and q^. They fulfill the transformation relations 






^7,7'' 

"7i7'' 



?;;7' = ?^,7', (53) 

where 7 and 7 are spin-orbitals with opposite spins. The coexistent superconducting- 
antiferromagnetic phase (SC+AF) can be defined in the following way 



Af^ = A^i^ = A? ^ 



.G 



-lA 



.G - aG 



A^^ = A^ ^ 



(54) 



Considerations analogical to those presented in subsection 2 lead to the conclusion that 
for both sublattices the non-diagonal variational parameters, Xfj, and Xfjr, that have 
to be used in the calculations, appropriate for the SC+AF phase, are the same as those 
listed in Table [2l This fact, and the degeneracy of our bands, allow us to apply f l35l) for 
both sets of renormalization factors (for A and B sublattices), as we have 



(55) 



where a represents the spin opposite to a. Now, we can write down the Hamiltonian 
Kga for the case of SC+AF phase 



5lo-,lcr 


— Q2a,2a 


= ^f. 


la — 


<ll 


2ct — 


Qa 


_A 

Q2a,la 


= 5lo-,2o- 


= ^t. 


la ~ 


t-a 


2a = 


Qa 


qLa = 


= (l2aA = 


(flcrB 


= Q2. 


tB - 


-^ql 





K, 



GA 



Yl <5ek(4aACk«aB + 4«<xBCk«aA) 



Ma 



Ml'a 



Qekuic^l^A^kl'aB + C^i^B^M'aA) 



- ^ Y^^qlfl^laA + qlfl^^laB) + 1^Y. ^thiKjJo + IT 5^ EIiAKm)0 



(56) 



kla 



h,h 



h,h 



where 

Q = Qm - Mi ■ (57) 

It should be noted that the sums in (156|) are taken over all L/2 independent k states. 
As before, we apply the SGA method which leads to the effective Hamiltonian with the 
statistical-consistency constraints of the form 



K.=K, 



GA 



Ac 



7^(y{fii,iaA - nuaB) - 2L5'| 



^2' 

klo 



OU 



-A_ 



-A. 



>- k 



iklA + ^k-lBj 
Ik-IA + ^klfi) 



LA 



A. 



k 

k 

^iq^nklaA + qlflMaB) - Lug 



LA" 



H.C. 
H.C. 



(5J 



I- k«o 
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Introducing now the eight-component composite operator 






I'^klcrA' '^k2o-A' '^klcrS' '^k2o-_B' C-klo-A, C-k2aA, C-klo-B, C-k2aB) , 



we can write down the effective Hamiltonian K\ in the following form 
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where the explicit form of the 8x8 matrix is 
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Diagonalization of ( 16T]) leads to the quasi-particle energies E'kifr (^ = 1? 2, 3, . 
corresponding grand potential function F\ per atomic site now has the form 
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(59) 
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8). The 



(63) 



As before, we minimize the Fx function to determine the values of the mean fields, the 
variational parameters and the Lagrange parameters. The necessary conditions for the 
minimum are again expressed by fHSl) and (H6|) . In the subsequent discussion we consider 
also the pure antiferromagnetic phase (AF), for which S*^,^ ^ but A+ = A_ = 0. The 
number of equations that need to be solved is different for different phases considered 
in this work. In Table III we show how many equations are included in (145|) and (Il6|) 
for all phases discussed. 
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Table 3. Number of equations that have to be solved in the case of ah considered here 
phases. To reduce the number of equations for particular phases we have used certain 
symmetry relations regarding the mean field parameters, the Lagrange multipliers, and 
the variational parameters. 



phase 


A 


Al+FM 


SC+AF 


AF 


NS 


FM 


num. of Eq. 


16 


17 


22 


12 


8 


13 



3. Results and discussion 



Equations (H5|) and ( H6|) have been solved numerically for all phases by means of the 
hybrdl subroutine from the MINPACK library, which performs a modification of the 
Powell hybrid method. The maximal estimated error of the procedure was set to 10~^. 
The derivatives in Eq. ( l45l) and (146|) were computed by using a 5-step stencil method 
with the step equal to x = 10~^. 

We concentrate now on the detailed numerical analysis of the phase diagram and the 
microscopic characterization of the stable phases. Having in mind that for 3d orbitals 
U' = U — 2J, one obtains the HF condition for the pairing to occur, U < 3J (see |22)). 
We discuss thus first and foremost the limit U < 3 J, as it allows for a direct comparison 
of SGA with the HF solution. In this manner we can single out explicitly the role of 
correlations in stabilizing the relevant phases. One should note that in the considered 
regime [U < 3J) we have a model with intraatomic interorbital attractions leading to 
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Figure 1. (Color online) Stable phases evolution vs. band filling. The 
superconducting gap parameter, magnetic moment and free energies as a function 
of band fining both for the HF and SGA, for J = 0.299: a, b, c and J = 0.4545: d, e, f. 
The results are for /3/j = 0.0. The shaded regions represent the stability regions of the 
respective phases according to the SGA calculations. In Figs, a and d we show only 
the free energies of stable phases. The arrows in a and d mark the transitions points 
between phases. 
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spin-triplet pairs. As the main attractive force is of intraatoniic nature, we focus here 
on the local (s-wave) type of pairing only. In other words, as we discuss the situation 
with no or small hybridization, the intersite part of the pairing can be disregarded. 

The calculations have been carried out for U' = U — 2J, U = 2.2 J, ksT/W = 10~^. 
This leaves us still with three independent microscopic parameters in our model: no, 
J, and Ph. All the energies have been normalized to the bare band- width W = 8\t\ 
(as we consider the square lattice with nearest neighbor hopping). For comparison, we 
also show the results calculated by means of the combined HF-BCS=HF approximation. 
This method is described in detail in our previous paper for the same model as considered 
here. We can also reproduce the HF results by using the Gutzwiller method described 
in this work and setting Xjj/ = 6jj'. 

In Fig. [T]we display the free energy, superconducting gaps, and magnetic moments 
for the two values J = 0.299 and J = 0.4545. As one can see from the free-energy plots 
(Figs. [1^ and[T]i), below some certain value of band filling, the pure superconducting 
phase of type A is stable for the SGA method. The increase of the number of electrons 
in the system, enhances the gap in this region (Figs. [Tb and[T^). Above the critical 
band filling Uc, the staggered moment structure is created and a division into two gap 
parameters (A+ and A_) appears, as can be seen in Figs. [TJ), [1^, Ht, and [If. In this 
regime the SC+AF phase becomes stable. 

When approaching half filling, both gaps gradually approach zero and for n = 2 
we are left with a pure AF phase, which is of Slater insulating type evolving towards 
the Mott-Hubbard insulating state with the increasing U. As the staggered magnetic 
moment is rising (with the increase of no), the renormalization factor is approaching 
unity (cf. Insets to Fig. [1^ and[T]i). This is a consequence of the fact that for large 
values of Sq, the configurations with two electrons of opposite spin, on the same orbital, 
are ruled out. 
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Figure 2. (Color online) The superconducting gaps (a) and the free energies (b) as 
a function of band filling for J = 0.4545 and /3^ = 0.1. The shaded regions represent 
stability of respective phases according to the SGA calculations. The vertical arrows 
mark the phase borders. 



Comparing Figs. [T^, [TJd, [Tt with Figs. [TJl, [T^, [If one sees that by increasing J 
we make the value of n^ smaller. However, the decrease in n^ is not as significant in 
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SGA as it is in the HF case. In general the results presented Figs. [TJo, [T]:, [T^, [Tf 
look similar from the qualitative point of view for both methods. For SGA, the onset 
of antiferromagnetically ordered phase appears closer to half filling than for the HF 
method. Another difference between HF and SGA is that for the former the staggered 
moment in the SC+AF phase is increased by the appearance of SC for the whole range 
of band fillings, whereas in SGA calculations the staggered moment is slightly stronger 
in the AF phase than in the SC+AF phase for a small region close to the half-filled 
situation (inset of Fig. [If). 
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Figure 3. (Color online) The superconducting gaps and magnetic moment as a 
function of the hybridization strength, /3/i, for n — 1.9, J — 0.4545, for the case of 
SC+AF and AF phases. 

Significant differences between HF and SGA can be seen in Figs. [It and [TJ'. 
While changing the band filling from to 2, in the case of SGA calculations we move 
consecutively through the regions of stability of NS (for J = 0.299), A, SC+AF phases, 
and for n = 2 we have pure antiferromagnetism. The situation is different in the 
HF approximation, where in between the regions of stability of A and SC+AF phase, 
we have also the stable Al+FM phase. It should be also noted that the free energy 
calculated in SGA is lower than the one for the HF situation, as one should expect, 
since the correlations are accounted more accurately in the former method. It is also 
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Figure 4. (Color online) The superconducting gaps in the SC+AF phase (a) and free 
energies of stable phases (b) as a function of Hund's coupling for n = 1.9 and 13^ = 0.1. 
The shaded region represent the stability of NS phase according to the SGA results. 
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Figure 5. (Color online) The superconducting gaps in the A phase as a function of 
the Hund's coupling for n = 1.0 (a-for [3h = 0.0 and c-for (3h = 0.1) and free energies 
of stable phases corresponding to SGA and HF approximations (b-for /3h = 0.0 and 
d-for /3ft, = 0.1). The shaded regions represent the stability of the NS phase according 
to the SGA. The vertical arrows mark the border points between respective phases. 
Insets: Bandwidth renormalization factor for /3/j = (upper) and (3^ =0.1 (lower). 



very interesting that having the system with U < 3J, no pure ferromagnetisni appears 
in this canonical model of itinerant magnetism. 

In Fig. [21 we present the results for the case with nonzero hybridization parameter, 
/3/i = 0.1. In this case there are no superconducting solutions below some certain value 
of the band filling (cf. Fig. |2^) and an extended region of NS stability occurs. The 
influence of the hybridization on the antiferromagnetically ordered phases is weak, as 
can be seen more clearly in Fig. |3l The changes in the superconducting gap and the 
magnetic moment triggered by the hybridization, are quite small even for larger values 
of/3/.. 

Next, we discuss the J dependence of the superconducting gap, the free energy and 
the magnetic moment for selected values of band filling. As in the case of n-dependences 
the gap parameters and the magnetic moments in both SGA and HF approximation are 
qualitatively similar. In Fig. IHwe can see that for n = 1.9 even the free-energy plots and 
regions of stability of certain phases are comparable for both calculation schemes used. 
For the quarter-filled case (cf. Fig. [5]) the Al+FM phase is stable above some value of 
J, according to the HF results. However, this is not the case in the SGA approximation, 
where the A phase has lower free energy even than the saturated ferromagnetic phase 
coexisting with superconductivity. Comparing Figs. [Sjo and [Sji (as well as [TJl and [JJo) 
one sees that the region of stability of the A phase narrows down in favor of the NS 
phase, due to the influence of hybridization. 
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It is important to check whether the itinerant magnetic phases are stable in the 
regime U' > J [U > 3J), i.e., when the superconductivity is absent in the HF 
approximation. For this purpose, in Fig. [6] we provide the band-filhng dependence 
of the free energy corresponding to stable phases for U = 4 J. Indeed, the paramagnetic 
and the magnetically ordered phases are stable for both methods of calculations. 
Therefore, for f/ > 3 J we recover the magnetic phase diagram for this model, which was 
considered originally only in the context of magnetism. The free energy of the saturated 
ferromagnetic phase calculated by the SGA is very close to the one obtained by the HF 
approach. This is again caused by the circumstance that in the saturated state all of 
the spins are parallel and the double occupancies on the same orbital are absent. In 
this situation, the intra-orbital Coulomb interaction is automatically switched off. It 
would be interesting to determine the stability of the coexistent phases in this regime 
{W > J). Work along this line is in progress. 
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Figure 6. (Color online) The free energies of the stable phases for SGA and HF 
methods for J ~ 0.4 and U = 1.6. The shaded regions in the inset mark the stability of 
certain phases according to the SGA approach. Note the appearance of ferromagnetic 
phase for L/ = 4J (i.e., for I] > 3J) in the filling range 1.45 -^ 1.75, sandwiched in 
between paramagnetic and antiferromagnetic phases. 



4. Conclusions and outlook 



4.I. Conclusions 

The principal purpose of this paper was to formulate a many-particle method which 
allows to investigate the spin-triplet real-space pairing in correlated system with an 
orbital degeneracy. To this end, we have carried out a detailed analysis using the 
statistically-consistent Gutzwiller approximation (SGA) for the two-band degenerate 
Hubbard model with the spin-triplet superconductivity and itinerant magnetism 
included, both treated on equal footing. The results were compared with those coming 
from the Hartree-Fock approximation amended with the Bardeen-Cooper-Schrieffer 
(BCS) approach. The obtained Hund's coupling and band filling dependences of the 
magnetic moment and the superconducting gap parameters are often similar from the 
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Table 4. Exemplary values of the order parameters, the chemical potential, the free 
energy, and the band renormalization factors corresponding to the considered phases, 
for two different sets of values of the microscopic parameters n and J. The underlined 
values correspond the stable phases. The numerical accuracy is on the level of the last 
digit specified. 
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n = 1.9 


parameter 


phase 


J = 0.299 


J = 0.299 
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0.1307664 
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0.1638992 


A_ 


SC+AF 
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s^ 
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0.6736373 
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0.8840340 


0.6421089 
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SC+AF 
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0.9211224 


Q 


AF 
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0.9293567 



qualitative point of view with those evaluated by means of the HF approximation. 
However, the stability regions of the considered phases are significantly different for 
the two applied methods. In SGA, the stable coexisting superconducting-ferromagnetic 
phase is absent while it appears in the HF approximation in a certain range of J and 
n values. Furthermore, the coexistence of the paired state with antiferromagnetism 
appears much closer to the half-filled situation in SGA than in HF approximation. For 
n = 2 the superconductivity disappears and only the pure antiferromagnetism survives; 
this state can be termed a correlated Slater-insulator state, which evolves gradually into 
the Mott-Hubbard insulating state with increasing U > 1 and Sq^ -^ 1/2. 

The influence of hybridization for both approximations is similar. With an increase 
of the Ph parameter, the region of stability of the superconducting type-A phase narrows 
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down in favor of the NS state. On the other hand, the antiferromagnetic phase is not 
affected in any significant manner by an increase of Ph. 

The band renorniahzation factors approach unity as the interaction constants J, 
U and U' tend to zero, what represents an additional test of our numerical results 
correctness. Generally, in the low-coupling limit our present results reduce to those 
obtained in HF approximation analysed by us in |22| . as it should be. 

It is important to emphasize that for both approaches the phase diagrams have 
been obtained for U < 3 J, i.e. for relatively low value of the Hubbard interaction U, or 
equivalently, for a relatively high value of the Hund's rule exchange integral. A complete 
analysis of the present model would require studying the stability of the spin-triplet 
superconductivity and its coexistence with magnetic ordering in the complementary 
regime U > 3 J, where the magnetism is favored against superconductivity. This regime 
has been the subject in a number of earlier papers [211 H5l Il6] . as then both the 
intraorbital, as well as the interorbital interaction is repulsive, and lead in a natural 
manner to magnetic ordering. 

4-2. Outlook: Extension and application to real systems 

In connection with the remarks provided above, we would like to characterize briefly 
the possibility of extending the present model ([T]) to the uranium systems in which 
superconductivity and ferromagnetism coexist in an unambiguous manner [36]. First 
of all, the magnetic moment in those systems, particularly in UGe2 and URhGe, is 
quite large, with an associated molecular field in the megagauss range, which most 
probably rules out any spin-singlet character of pairing (note that the Curie temperature 
(Tq) to the superconducting transition temperature (Ts) ratio reaches in the uranium 
compounds the value Tq/Ts ~ 10^). In spite of those circumstances, our solution 
does not provide any extended regime (for the studied parameter range) for the 
ferromagnetism-spin-triplet superconductivity coexistence. Instead, in a wide range 
of band fillings, the coexisting SC+AF phase is stable (cf. Fig. ^ and[I]F), as well as 
the pure spin-triplet superconducting phase of type A (the equal-spin-paired phase). 
The pure A phase seems to be realized in Sr2Ru04, though then a detailed three-orbital 
structure of the order parameter seems to be relevant [lO]. A direct application of 
our SGA scheme to a realistic three-band system is more involved, as the number of 
parameters to minimize would lead to a computing time-consuming procedure, but still 
possible to tackle. 

The extension of the present model to the uranium system such as UGe2 would 
require considering orbitally degenerate and correlated 5/^ — 5/'^ quasi-atomic states 
due to U and hybridized with the uncorrelated conduction band states. This means that 
we must have minimally a three-orbital system with two partially occupied 5/ quasi- 
atomic states (so the Hund's rule becomes operative) and at least one extra conduction 
band. Such situation may lead to a partial Mott-localization phenomenon, i.e., to a 
spontaneous decomposition of 5/"^ {n > 1) configuration of electrons into the localized 
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and the itinerant parts [H]. In such a situation, it is possible that the locahzed electrons 
are the source of ferroniagnetisni, whereas the itinerant particles are paired P . This is 
not the type of coexistent phase we have in mind here, since in the model considered by 
us all the system electrons are indistinguishable in the quantum-mechanical sense. 

These considerations lead to the conclusion that one would require minimally a 
periodic Anderson model with degenerate 5/" states, to mimic the uranium-based 
ferromagnetic superconductors. This variant of the multiple-band model is also very 
useful in the discussion of heavy-fermion compounds. Moreover, in the systems 
represented by this model, the coexistence of antiferromagnetism and superconductivity 
has been shown to appear in both experiment |18] and theory |19]. More specifically 
for the systems UPts and UNi2Al3 the coexistence of the spin-triplet pairing and 
the antiferromagnetism has already been suggested to appear |50| IST] . although not 
elaborated in any detail. One specific feature should be mentioned. Namely, in the 
situation when we have antiferromagnetic superconductor, then there is also a strong 
theoretical indication that there is a spin-triplet component even for the pure spin-singlet 
mechanism of pairing |27tH9]. The spin-triplet important admixture results simply from 
a decomposition of the system into two sublattices with staggered magnetic moment. 
These and related features must be taken into quantitative analysis before any realistic 
consideration of concrete systems is carried out. 

In connection with the whole discussion, it is intriguing to ask if a symmetric model 
system of the type exemplified by Hamiltonian ([1]) could be experimentally realized in 
the optical lattice. Some model systems (e.g., the Hubbard model system) have been 
experimentally achieved in this manner |42| . 

In relation to the spin-triplet real-space pairing induced by the Hund's rule, one 
should also mention the spin fiuctuations (SF) as a possible mechanism of spin-triplet 
pairing in both magnetic |l2] and liquid '^He systems |l3]. Within the present approach 
the spin fiuctuations should be treated as quantum fiuctuations around the present self- 
consistently renormalized mean field state [Hj. The real-space and the spin-fiuctuation 
contributions may become of comparable magnitude in the close vicinity of the quantum 
critical point, where the ferro- or antiferro- states disappear under e.g. pressure. This 
is, however, a completely separate topic of studies. 
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